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ABSTRACT 



Context. A positive power-law trend between the local surface densities of molecular gas, S ?aj , and young stellar objects, £*, in 
molecular clouds of the Solar Neighbourhood has recently been identified by Gutermuth et al. How it relates to the properties of 
embedded clusters, in particular to the recently established radius-density relation, has so far not been investigated. 
Aims. We model the development of the stellar component of molecular clumps as a function of time and initial local volume density. 
Our study provides a coherent framework able to explain both the molecular-cloud and embedded-cluster relations quoted above. 
Methods. We associate the observed volume density gradient of molecular clumps to a density-dependent free-fall time. The molecular 
clump star formation history is obtained by applying a constant star formation efficiency per free-fall time, e/f. 
Results. For volume density profiles typical of observed molecular clumps (i.e. power-law slope m —1.7), our model gives a star-gas 
surface-density relation of the form oc S 2 aj . This is in very good agreement with the relation observed by Gutermuth et al. Taking 
the case of a molecular clump of mass Mq = 10 4 M Q and radius R a 6pc experiencing star formation during 2Myr, we derive what 
star formation efficiency per free-fall time matches best the normalizations of the observed and predicted Eg flJ ) relations. We find 
ejf ca 0. 1 . We show that the observed growth of embedded clusters, embodied by their radius-density relation, corresponds to a surface 
density threshold being applied to developing star-forming regions. The consequences of our model in terms of cluster survivability 
after residual star-forming gas expulsion are that due to the locally high star formation efficiency in the inner part of star-forming 
regions, global star formation efficiency as low as 10% can lead to the formation of bound gas-free star clusters. 

Key words, stars: formation — galaxies: star clusters: general — ISM: clouds — stars: kinematics and dynamics 
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1. Introduction 

Properties of star-cluster forming regions are crucial in deter- 
mining how the nascent cluster dynamically responds when the 
gas left unprocessed by star formation is driven out due to stel- 
lar feedback. Does the cluster survive as a bound entity - albeit 
depleted of a fraction of its stars, or is it fully dispersed into 
the field? Pivotal properties influencing cluster-survival likeli- 
hood include the star formation efficiency of cluster-forming re- 
gions and their potential well depth. The latter can be quantified 
through their mass-radius relation, which several authors report 
to scale as m ec i oc r 2 t , where m ec i and r ec i are the embedded- 
cluster stellar mass and radius, respectively (lLada & LadaLl200"3t 
lAdams et al 1 I2006I) . We stress that this relation refers to the stars 
only, i.e. m ec i does not account for the unprocessed gas. 

This mass-radius relation is one of constant mean surface 
density and can stem from the method applied to measure the 
embedded-clus ter radius and the stellar mass it contains. As 
pointed out by lAllen et all J2007), if the mass and radius of a 
cluster are defined based on a surface density threshold, they are 
necessarily sensitive to the adopted density cut-off. This also de- 
termines what mass fraction of the stars is hosted by the 'halo' 
surrounding the star cluster, as opposed to the cluster itself. 
For a given star-forming region, the higher the surface density 
threshold, the smaller the cluster mass and radius, the larger the 
surrounding 'halo'. This renders an accurate definition of the 
cluster-forming region problematic. 



Besides, it is worth stressing that the dynamical response 
of a cluster to the residual star-forming gas expulsion depends 
on the cluster-forming region properties at the onset of gas 
expulsion. This prompts another key question when quantify- 
ing cluster-forming region properties: is an observed embedded 
cluster caught in the process of turning its gas into stars, or has it 
just reached the end point of star formation with the in tra-cluster 
gas about to be expelled? iHiguchi et al.l d2009l l20loh illustrate 
the problem well. Their C ls O mapping of 14 molecular clumps, 
each forming a star cluster, shows a sequence of star formation 
efficiencies. The highest efficiencies are found for clusters as- 
sociated to C I8 Oemission holes, highlighting that gas dispersal 
has started therd3 Therefore, an observed star formation effi- 
ciency is not enough to assess how the corresponding cluster 
will respond to gas expulsion. In Section l4.1.2l we will suggest 
an alternative explanation to gas dispersal to explain the C 18 0- 
emission holes of some cluster-forming clumps. 

The characterization of cluster-forming regions has made a 
significant leap forward thanks to the Spitzer Space Telescope. 
Spitzer has probed a wide spectrum of stellar densities in star- 
forming environments of the Solar Neighbourhood - from rel- 
ative isolation to star clusters (although the Trapezium region 
of the Orion Nebula Cluster remains unresolved with that facil- 
ity). 5/7/fzer-surveys have led to a new picture, that is, star clus- 



1 An alternative explanation to the C 18 depletion, however, could 
be that C 18 is frozen out on grain surfaces. Observations of the less- 
affected nitro gen-containing spec ies NH3 and N2H + should help clarify 
the question dBergin et alll2002h 
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ters are the emerging stellar peaks of wider star-forming region s 
(lAllen et all 120071: Evans et all 120091: IGutermuth etall l201ll) . 
Understanding the physics of cluster formation and the prop- 
erties of star-forming regions at large are therefore now tightly 
entwinned topics. 

IGutermuth et alJ d201 ll) identifies a positive power-law trend 
between the local surface densities of molecular gas and young 
stellar objects (YSOs) in eight molecular clouds located within 
1 kpc from the Sun. On the average, the relation between the lo- 
cal YSO surface density, "Lyso, an d the local gas surface den- 
sity, Y, g , follows 1.yso - 10~ 3 Sg, with a - 2 and the surface 
densities in units of M Q ■ pc~ 2 . Their result is based on map- 
ping Sp/fzer-identified YSO spatial distributions with a n ,h near- 
est neighbour scheme, and tracing the molecular gas by near- 
infrared extinction. In most molecular clouds of the sample, the 
po wer-law scaling i s affect ed by a large scatter (~ 1 dex, see fig. 9 
in IGutermuth et all 1201 ll) . It is tightest in the Ophiuchus and 
MonR2 clouds where the power-law index a is 1.9 and 2.7, re- 
spectively. One of the purposes of this contribution is to show 
how this star formation law and the physics of cluster formation 
are related. 

There is increasing evidence fo r a bimodal cluster forma- 
tion process. iDa Costae t al. (2009) find that the half-light ra- 
dius distri bution of extragalactic ol d globular clusters is bimodal. 
Similarly, iBaumgardt et al.l (1201 Ol) identify two populations of 
globular clusters in the outer halo of our Galaxy based on the ra- 
tio of their half-mass to Jacobi radii. They suggest this bimodal- 
ity to be an imprint of cluster formation, rather than a conse- 
quence o f 13Gyr of dynam ical evolution in the Galactic halo 
(see also lElmegreenl 120081) . In ad dition, a bimoda l cluster for- 
mation agrees with the finding of iPfalznerl (120091) that, in the 
Galactic disc, star clusters younger than 20Myr unfold along 
two distinct sequences in the space of cluster volume density 
versus cluster radi us (her fig. 2). She co ins them the starburst 
and leaky clusters. llMafz Apellanizl d2001l) identifies a similar di- 
chotomy in the st ructural properti es of a sample of nearby extra- 
galactic cl usters. 1 Pfalzne r (2009) also notes that the embeddded 
clusters of lLada & Ladal (120031) define a precursor sequence to 



the leaky clusters (her fig. 4). That is, these embedded clusters 
seem to be the leaky-cluster progenitors, with starburst cluster 
precursors still to be identified. 

In a follow-up study, Pfa lzner! d201 ll) shows that the 
embedded-cluster sequence obeys p ec \ oc r~\ , where p ec i is the 
cluster mean volume density (her fig. 2). This is reminiscent of 
a relation of c onstant mean surfac e de nsity, i.e. PeA oc r~ \, as 
put forward by lAdams et all d2006l) and lAllen etall d2007f) ! She 
presents the embedded-cluster sequence as a growth sequence, 
that is, clusters are still in the process of building their stellar 
mass and their properties may therefore not be representative 
of the conditions at gas expulsion onset. The embedded-cluster 
scaling p ec \ oc r~V 3 equates with m ec i oc r 1 ']. This implies that, 
as the embedded-cluster (stellar) mass m ec / increases with time, 
so does the radius r ec i, while the mean volume density p ec i de- 
creases. In other words, in this scenario, star formation propa- 
gates outwardly. The cluster central regions form first, and outer 
shells of stars are added with time. This is an extremely ap- 
pealing concept for the following reason. The timescale rele- 
vant for star formation is the free-fall time of the s tar-forming 
gas, T ff dElmegreenl [20071: iKrumholz & TanL 120071) . Since r// 

— 1/2 

scales with the gas volume density as Tff oc p and since 
molecular clumps have radial density gr adients dBeuther et all 
l200l iMuller etall l200l IPirogovL [2009h . the gas free-fall time 
is shorter in the denser inner regions of molecular clumps than 



in their outskirts. That is, star formation is faster closer to the 
clump centre than towards the clump edge. 

Here comes a subtle difference, however. Is star formation in 
the outer regions of a molecular clump delayed, or does star for- 
mation start at the same time all through the molecular clump, 
albeit at a slower rate at la rger distance fr om the clump cen- 
tre? The scenario devised bv IPfalznerl d201 ll) clearly fits the first 
hypothesis. In contrast, a molecular clump forming stars on a 
radially-varying timescale would match the second hypothesis. 

Considerable efforts have been dedicated to the hydrody- 
namical simulations of star formation in molecular gas (e.g . 
[ Klessen. Burkert & Batelll998l:lBate. Bonnell & BrommL 120031: 
Bonnell. Clark & Batd I2008T These simulations are compu- 
tationally expensive and hence limited in terms of the mod - 
elled gas mass (e.g. 10 4 M Q , iBonnell. Clark & BateL 120081) . 
Parametric semi-analytical studies remain therefore useful, es- 
pecially to browse the parameter space extensively. 

In this contribution, we build on the concept of star for- 
mation efficiency per free -fall time originally introduced by 
IKrumholz & McKea d2005l) . namely, the fraction of an object's 
gaseous mass that is processed into stars over one free-fall time 
at the mean density of the objecfl- We generalize their approach 
by defining a local free-fall time. That is, we present a new 
model for star-forming regions which hinges on the radially- 
varying free-fall time of molecular clumps. We show how it pro- 
vides an elegant and coherent picture accounting for both the 
local B star formation l aw Zyso x ^i a s recenly highlighted by 
IGutermuth et al.l d2011[). and the grow th sequence of embedded 
clusters postulated bv IPfalznerl d201 ll) . 

Our paper is organised as follows. In Section [2] we present 
the model and show how it e xplain s straightforwardly the scal- 
ing law of IGutermuth et al.l d201ll) . In Section [3] we fold the 
model with a surface density threshold. This equates with ac- 
counting for a contaminating stellar background against which 
the star-forming region is projected and observed. We deduce 
the corresponding r ec i-p ec i relation, where r ec i and p ec \ are the 
radius and mean volume density of the part of the star-forming 
region which 'emerges above' the stellar background. We find 
the predicted r ec [-p ec \ relation to be in good ag reement with the 
embedded-cluster sequence of IPfalznerl (1201 ll) . Section [4] gives 
an overview of the model consequences in the framework of 
cluster gas expulsion. Finally, we outline some future-work di- 
rections in Section|5]and present our conclusions in Section|6] 



2. Free-fall time driven star formation 

We start by studying how a spherical gas clump with a radial 
density gradient builds its stellar component as a function of 
time and space, under the assumption of a constant star forma- 
tion efficiency per free-fall time, ejf. Because of the density gra- 
dient, the free-fall time must be defined locally, i.e. Tff(r), and 
is an increasing function of the distance r from the clump cen- 
tre. Consequently, star formation proceeds more quickly in the 
clump central regions than in its outskirts. This leads to a density 
profile for the stellar component steeper than the initial density 
profile of the clump. 



2 Note that lKrumholz & McKed J2005T) use the terminology 'star for- 
mation rate per free-fall time', or SFRyy for short. We prefer to refer to 
the 'star formation efficiency per free-fall time' and denote it ejf. 

3 We refer to the Gutermuth et al. ( 20 11 ) star formation la w as a local 
law to make it distinct from the global star formation law of Kennicutll 

1 which builds on galaxy-integrated surface densities. 
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Prior to going into detailed numerical simulations, we 
present a simple analytical approximation which allows to grasp 
the relevant physics. 

2.1. Analytical Insights 

Assuming spherical symmetry, the volume density profile, po(r), 
of a molecular clump of mass Mo, radius Ro, and density index 
po, obeys: 

3-po M 



po(r)=k p , r- po = 



-Po 



An 



(1) 



where r is the distance from the clump centre. The subscript 
'0' refers to the gas properties prior to star formation (i.e. at 
t — 0: the clump is made of gas only). The factor k p o comes 
from integrating the density profile over the clump volume, i.e. 

Mo = j Q R ° 4nr 2 po(r)dr. The density index of star-forming molec- 
ular clumps is observed to ra nge from = 1.5 to - 2.0. We there- 
fore a ssume 1.5 < po < 2.0 dBeufher et al.L l2002t iMuller et all 
I2001 FI. 

In this model, we assume that star-forming molecular 
clumps experience neither significant outflows or inflows, nor 
do they contract or expand during star formation. We also 
neglect the potential migration of YSOs after their forma- 
tion. As a result of these assumptions, at any time t after the 
onset of star formation, the clump total mass (gas + stars) 
is preserved, i.e. M*(f) + M g (t) = Mo, as well as the clump 
radius Ro and the spatial distribution of the clump mass, i.e. 
pi,(t,r) + p g (t,r) = po{r). The subscripts 'g' and refer to 
the properties of, respectively, the unprocessed gas and stellar 
component at any time t > 0. As an example, we write p g and 
po for the current and initial gas volume densities. Finally, we 
consider that the star formation efficiency per free-fall time, e/f, 
is constant. In this class of model, star formation takes place 
all over the clump - albeit at a slower rate in the outskirts than 
in the centre. The limiting radii of the clump and of its stellar 
component are therefore equal: Ro — R*. In what follows, we 
shall refer both radii as R. 

Now let us consider that, at any time t, a fraction etf of the 
gas mass at radius r is turned into stars every local free-fall time. 
Because of the volume density gradient of the molecular clump, 
star formation proceeds more quickly in its centre than in its 
outskirts. As a result, the volume density profile of the stellar 
component built by the clump, p*(f, r), is steeper than po(r). We 
insist that this stems from the shorter free-fall time at smaller 
radius, not from a higher e/y in the clump centre. The free-fall 
time of the gas at radius r and time t obeys: 



T ff (t,r): 



3n 



32Gpg(t,r) 



(2) 



with G the gravitational constant and p g (t , r) the volume density 
profile of the unprocessed gas at time t. 

Let us assume that the volume density profile of the stellar 
component is a power law too, with a constant density index q: 

3-?M*(f) _ a 



An R^-D 



(3) 



4 St rictly speaking, th e dens ity indices inferred by iBeuther et alj 
(2002) and Miil ler et alj l l2002j) characterize star-forming molecular 
clumps. They may thus be different from po which, in this model, is 
the gas density index prior to star formation. We shall come back to 
this point in Section l531 
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Fig. 1. Top panel: Initial gas volume density profiles (filled sym- 
bols) and their associated steeper stellar density profiles (open 
symbols). The latter are obtained from the numerical model one 
million years after star formation onset, i.e. t — 1 Myr, for three 
distinct po density indices (see key). The spherical molecular 
clump has a mass Mo - 10 4 M Q , a radius R =2 6 pc, and the star 
formation efficiency per free-fall time is e/f = 0.1. Middle panel: 
Stellar density profiles from the numerical model (lines with 
open symbols) for the times t and density indices po quoted in 
the key. Each symbol-free line is the corresponding analytical 
upper limit (Eq. [8}. Bottom panel: Time evolution of the gas 
volume density profile for po = 1.7. The star-formation-driven 
depletion of the gas in the central region of the molecular clump 
- where the free-fall time is the shortest - is clearly highlighted. 
The symbol-free dotted line over the range 0.2-1 pc has a slope 
of -1.3, i.e. shallower than the initial slope —1.7. 



As we shall see from the numerical modelling, this is a realistic 
approximation (see top and middle panels of Fig.[T]and top panel 
of Fig. |5j. In Eq. [3] M*(f) is the total stellar mass contained by 
the clump at time t. As previously, the factor k Pt + stems from 
integrating the density profile over the entire clump volume, i.e. 

M *( f ) = f R 47Tr 2 p*(t,r)dr. 
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Because the stellar component is steeper than the initial gas 
density profile, we have q > po- How much steeper is q com- 
pared to pol Let us assume for a moment that, at any time t 
and radius r, the gas mass keeps dominating the stellar mass. 
That is, in p+(t,r)+p g (t,r) = po(r), we assume p±{t,r) «p g (t,r), 
which leads to p g (t, r) = po(r), i.e. the gas density profile does not 
evolve significantly. Therefore, the local free-fall time (Eq.|2]) is 
approximately constant too: 



Tff(t,r) =* T ff (t = 0,r) . 



(4) 



Now consider a shell of thickness dr at a distance r from the 
clump centre and of initial gas mass dmo(r) = 4nr 2 po(r)dr. 
Under the assumption of a constant free-fall time (Eq. @}, the 
stellar mass formed by this shell at time t, dm+(t,r), follows 
from: 

dm+it, r) < e f f — —-dmo(r), (5) 

T ff (t = 0,r) 

where f/T//(? = 0, r) is the number of free-fall times elapsed 
since star formation started. The rhs of Eq. [5] defines an upper 
limit to the actual mass in YSOs in the shell at time t, for two 
reasons. Firstly, as time goes by, the star formation efficiency per 
free-fall time is applied to an ever lower gas mass, i.e. dm g (t, r) < 
dmo(r). Secondly, the steady decrease of the gas mass lengthens 
the local free-fall time, i.e. we have Tff(t,r) > Tff(t = 0,r) in- 
stead of Eq.|4] Equation [5] can be rewritten: 



P*(t,r)<e ff - 



T ff (t = 0,r) 
Combining Eqs.[2]and[6]leads to: 



Po(r), 



32G 

\/-^r boOO] 



3/2 



Introducing Eq.Q]then gives: 



P*(t,r)< 



32G 



t-k 



3/2 . r -3p /2 



(6) 



(7) 



(8) 



The rhs of Eq.[8]is shown as the symbol-free lines in the mid- 
dle panel of Fig.Q] The lines with open symbols show the exact 
solutions obtained either numerically (Section 12.2b or analyti- 
cally (Section [2.3t and which account for the time variations of 
the gas mass and of the local free-fall time. The panel highlights 
clearly that the analytical approximation provides an upper limit 
to the actual stellar density profile. We assume Mo =* 1O 4 M , 
R^dpc and e// = 0.1 (values discussed in Section |2~2| ). Star- 
formation durations t and density indices po are given in the key. 

The rhs of Eq. [8] works best as an analytical approximation 
for p*(f,r) in the clump outskirts where the local free-fall time 
is the longest. That is, the hypothesis p g (t,r) =* po(r) remains 
valid for the time-span of our simulations (up to 2.5 Myr) 
and Eq. [8] can be seen as an equality (rather than an upper 
limit). In the inner regions, however, the shorter free-fall time 
implies a gas depletion quicker than in the outskirts. This 
eventually results in p g {t,r) < po{r), > T//(? = 0,r) (star 

formation slows down), and the rhs of Eq. [8] gives a firm upper 
limit. This differential behaviour between the inner and outer 
regions means that the actual stellar density profile is shallower 
than given by the rhs of Eq. [8] (compare the symbol-free 
lines and the lines with open symbols in the middle panel of 
Fig. Q] especially for the po = 1.9 model). In other words, the 
stellar density profile is steeper than the initial gas density 
profile (Eq. [7j by at most a factor 1.5, e.g. an initial gas density 



index po = 1.7 gives rise to a star density index q < 3po/2 = 2.55. 

We now have to convert the volume density profiles into 
(projected) surface density profiles, since observer-retrieved 
quantities are surface densities. For power-law volume density 
profiles, the surface density profiles are shallower than their vol- 
ume counterparts by 1 dex. Therefore, the surface density pro- 
files of the gas, initially and at time f, and of the stars obey: 



2o(s) = fc,o s 



-po+1 = 3-/7Q Mp __ pn+1 
27T Rl-Po 



E ? (M) = ^-' +1 = ^^-' +1 



2n B?-p 



and 



£*(/,5)=&: 



3-qMM 
2n /?3-« 



(9) 



(10) 



(ID 



respectively. The factors £x,(0,g,*) come from integrating the 
surface densities over the whole molecular clump surface: 

rR 

M{o, g ,+) — J T.(o, g ,-k)(s)-2ns-ds, where s is the projected distance 
from the clump centre. That is, s is a two-dimensional distance 
on the plane of the sky, while r is a three-dimensional distance. 

With the surface density profiles (Eqs.l9lfT0land[TT1i derived 
above and an estimate of the star density index (q < 3po/2), we 
are now ready to infer the local star formation law Z+(f, s) oc 
Y, g (t,s) a predicted by the analytical approximation. Eliminating 
the two-dimensional distance s between Eqs.[lO]and[TT] we ob- 
tain: 



(3-q)MM 



■{2nR 2 )"' [ ■(ZgCf.s))^ . (12) 



[(3 



Under our assumption that the gas mass dominates the stellar 
mass all through the clump, Eo(^) = %g(t, s) + ^*(t, s) - Z g (t,s) 
and we substitute £o(s) to s) in Eq.Q/2] This leads to: 



(3- 9 )M*(r) 



[(3 



■po)MoV'°- { 



■{2nR / -)Pa 



^\ Df\-V 



(13) 



Considering the index a of this relation 



p - 1 po - 1 ' 



(14) 



it immediately appears that a density profile steeper for the stel- 
lar compone nt than for the gas (i.e. q > po) is conducive to a > 1, 
as found by iGutermuth et al.l (1201 ll) . Were the density indices 
for the gas and stars identical (i.e. q = po, equivalent to no ra- 
dial variations of the local star formation efficiency; see also 
Section 14.1.2b , a would be unity. There is therefore a causal 
link between (i) the slope of the relation between the gas and 
YSO surface densities (Eq. [T3l . and ( ii) how much steeper the 
volume density profile of the stars is compared to that of the 
gas (Eq. |7]i. Adopting po - 1.7 and q < 3po/2 = 2.55 leads to 
a < 2.2, i n fair agreement with th e mean power-law trend ob- 
served by IGutermuth et aD (1201 ll) . This quantitative result pro- 
vides us with a strong incentive to refine the (L* , S ? ) relation by 
means of numerical simulations. 
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2.2. Numerical Model 

In the numerical model, we discretize the spherical molecular 
clump into successive shells defined by their three-dimensional 
radius r, local volume density p g (t,r) and local free-fall time 
T/f(t,r) at time t. Every instantaneous local free-fall time, a frac- 
tion £ff of the gas mass is removed and added to the stellar con- 
tent: 

(ti-ti-i) 

dmJti,r) = dmJti-\,r) - eff ■ dmJti-\,r). (15) 

Tff(t t -i,r) 

(ti-ti-i) 

dm+(ti, r) = dm+(ti-i, r) + e/f • dm g (ti-u r). (16) 

The local gas free-fall time (Eq. 13 is then updated. As the gas 
gets depleted, star formation slows down, an effect unaccounted 
for in the analytical approximation of Section [2~Tl 

The numerical results are shown in FigQ] for the case of 
a spherical gas clump of mass Mq = 1O 4 M , radius R — 6pc 
and star formation efficiency per free-fall time e/f = 0.1. Our 
choice of the mass and radius stems from inspecting the extinc- 
tion map of the MonR2 molecula r cloud and its spatial d istri- 
bution of 5p/fzer-identified YSOs ( iGutermufh et al.L 1201 ll their 
fig. 1). We focus on the largest concentration of YSOs (at a right 
ascension of 06 A 08 m and declination of =* -6°24') for which 
we adopt a radius R — 6 pc so as to include the lowest YSO sur- 
face densities. We adopt Mq = 1O 4 M as a rough estimate of the 
initial gas mass within R = 6 pc given that this region represents 
a significant fraction of the MonR2 cloud wh ose total mass is 
25,8OOM (table 1 in iGutermuth et all I20TH) . We plan to map 
our model onto detailed observational data in the near future. 
In the mean time, we stress that the estimate of e/f we adopt 
to make our model match the averaged local star formation law, 

= 10~ 3 S^, depends on the chosen values for Mq and R (see 
below and top panel of Fig. [3]). 

As for the upper limi t on the time t elapsed si nce star forma- 
tion onset, we note that IGutermuth et al.l d201 ll) study encom- 
passes mostly Class I protostars and Class II pre-main-sequence 
stars. Given that the respective average lifetimes of the Class 
I and Class II phases are 0.5 Myr and 2 Myr dEvans et all 
120091) . we adopt 2.5 Myr as the time-span of our simulations. 

The top panel of Figfl] illustrates the initial volume density 
profile of the gas for different density indices: po = 1.5,1.7 and 
1.9 (plain symbols: red squares, green circles and blue triangles, 
respectively). Also depicted are the volume density profiles of 
the built-in stellar component at t — 1 Myr (open symbols with 
identical symbol/colour-coding), thereby highlighting the steep- 
ening of the stellar density profile compared to the initial spatial 
distribution of the gas. The gas-to-star steepening equates with a 
local star formation efficiency higher in the clump central region 
than in its outskirts. We come back to this point in Section !?. 1.21 
The middle panel shows the stellar density profiles for the gas 
density indices po and times t quoted in the key. Each profile is 
shown along with its upper limit predicted by Eq. [8] The differ- 
ence between the analytical approximation (symbol-free lines) 
and the numerical model (open symbols) is stronger when the 
free-fall time is short, e.g. in the clump inner regions, especially 
for a steep density profile (high po). The difference also gets 
higher for longer star-formation durations t. Nevertheless, the 
comparison between both demonstrates the excellence of Eq. [8] 
in providing a back-of-the-envelope estimate of the stellar den- 
sity profile. The bottom panel of FigQ] depicts the time evolu- 
tion of the gas density profile when po = 1.7. Its flattening in 
the clump centre contrasts markedly with the absence of evo- 
lution at the clump edge where the initial free-fall time is the 



longest. The instantaneous gas density index, p, thus becomes 
smaller/shallower as t increases. 

2.3. Analytical solution 

The actual time-evolution of the star and gas density profiles 
can also be obtained analytically. Equations [T31fT6l correspond 
to separable first order differential equations. Using Eq.|2] they 
can be rewritten as : 



dp g (t,r) 
dt 



£ ff 



T ff (t,r) 



32G 



Pg( t ' r ) = -\l-j—- e ffPg( t > r ) 



3/2 



dt T ff (t,r) 
and their solutions are: 



Pg(t,r) 



32G 
~3n 



Pg(t,r) 



P0(rT U2 + 



■e f fp g (t,r) 



3/2 



P*0,r)=p (r)- 




(17) 



(18) 



(19) 



(20) 



where the clump density profile po(r) is given by Eq. Q] 
Equations [19] and [20] respectively correspond to the lines with 
plain symbols in bottom panel of Fig. Q] and to the lines with 
open symbols in top and middle panels of Fig. Q] obtained 
numerically above. 

Before going any further, we remind that our model builds 
on the local free-fall time, that is, it quantifies the local (at the 
radius r) collapse of the gas. The global collapse of the star- 
forming molecular clump (i.e. its collapse to its centre) may also 
be relevant to its evolution. The free-fall time to the clump cen- 
tre depends on the mean volume density enclosed within the ra- 
dius r, po(< r), and may therefore be shorter than the local free- 
fall time, especially for steep density pro files (e.g. the outskirt s 
of some molecular clumps; see fig. 2 in iBeuther et all 120021) . 
That is, were the clump evolution driven by gravity only, the gas 
would collapse to the clump centre faster than it does locally. 
The significance of the global collapse with respect to the lo- 
cal one depends therefore on how much supported the gas is by 
the combined effects of turbulence, magnetic fields, ... A study 
of that aspect is beyond the scope of the present paper and we 
assume for now that the star-formation history of the molecular 
clump is solely driven by the local collapse of its gas. As already 
quoted in Section 12.11 once the YSOs have formed at radius r 
through the gas local collapse, we neglect their potential migra- 
tion towards other regions of the molecular clump. The actual 
behaviour of newly-formed stars depends on their velocities at 
formation. If their initial velocity is sub-virial, once decoupled 
from the gas, they will fall to the cluster ce ntre where they will 
revirialize (see e.g. lGirichidis et all l2012ab . 

2.4. Surface densities and the local star formation law 

Following the analytical aproximations obtained in Section 12.11 
(Eqs [TOlfTTI) . we now derive the numerical solutions for the gas 
and star surface density profiles. To derive the exact surfa ce den- 
sity profiles, we build on eq. 6 of iParmentier et all d201 ll) which 
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Fig. 2. Surface density profiles from the numerical model. Lines 
with circles from top to bottom: the molecular clump gas ini- 
tially, the unprocessed gas and the stellar component at t — 1 Myr 
(see key). The molecular clump is the same as in Fig.Q] and the 
initial gas density index is po = 1.7. The solid symbol-free line 
depicts a slope of -0.7, that is, the slope expected for the initial 
surface density profile of the gas when po = 1.7. We note that 
in the clump outskirts, the actual surface density profile is be- 
coming increasingly steeper than the power-law approximation. 



integrates a power-law volume density profile into the projected 
mass enclosed within an aperture of radius s ap (see also their 
fig. 3). This provides straightforwardly the projected mass within 
a circular corona of radius s ap and thickness ds ap . Gas and star 
surface density profiles are shown in Fig.[2]for the same molec- 
ular clump as previously, po = 1.7 and the times t quoted in 
the key. The local star-formation law predicted by our numer- 
ical model can now be derived and compared with the scaling 
observed bv lGutermuth et all d2011l) . 

Figure [3] depicts the local surface density of YSOs, 
2*, in dependence of the local surface density of the ob- 
served/unprocessed gas, Hg. Both surface densities are in units of 
M • pc~ 2 . The (b lack) dotted line in the to p and middle panels is 
the scaling law of lGutermuth et all d201 ll) . i.e. = lO -5 ! 2 . The 
top panel shows the numerical model for po = 1 .7 at t — 0.5 Myr 
and t = 2.0 Myr (lines with open and plain symbols, respec- 
tively). The agreement between the model at t = 2.0 Myr and the 
mean observed scaling is excellent, both in terms of slope and 
normalization. This agreement is parameter-dependent, how- 
ever, as a molecular clump more massive (hence denser), or 
a longer time-span, or a higher star formation efficiency per 
free-fall time, all would shift the normalization upwards. The 
t = 2 Myr model bends upwards at high surface density. This is 
due to the flattening of the gas density profile in the clump cen- 
tre, that is, p decreases which in turn increases the index a of 
the local star formation law (Eq. [T4J, see also bottom panel of 
Fig.0. 

It is worth keeping in mind that = 10 3 £ 2 defines an 
average relation. The YSO-vs-gas relation observed for each 
molecular cloud exhibits a conside rable scatter (see the panels 
of fig. 9 in iGutermuth et all 1201 ll) . This is because a molecu- 
lar cloud consists of several molecular clumps each with its own 
mass, radius, density index and duration of the star-formation 
process. Therefore, a molecular cloud combines several corre- 
lations, each with its own normalization, corresponding to the 
several molecular clumps it contains. This leads to a trend in 
the observed (E ? ,£*) space, rather than the one-to-one rela- 



tion predicted for a single molecular clump and depicted in 
Fig. [3] Other reasons for the spread in the observed (£ ? ,Z*) re- 
lation include: Y SO migration and gas d ispersal via stellar feed- 
back process es dGutermufh et al.|, 201 ll) . non-spherical molec- 
ular clumps. IGutermuth et al.l (1201 ll) note the scatter to be the 
smallest in the MonR2 and Ophiuchus molecular clouds, which 
they represent as the green and blue parallelograms in their fig. 9. 
They are reproduced as the orange and brown parallelograms in 
the top panel of our Fig. [3] 

The middle panel of Fig. [3] illustrates how the model re- 
sponds to varying the gas density index po- For steeper profiles, 
the clump inner regions become denser at the expense of the out- 
skirts. This stretches the model towards both lower and higher 
surface densities, although the normalization is hardly affected. 

The bottom panel completes the top panel with the time evo- 
lution of the gas and star surface densities at four projected 
distances s from the clump centre (black dotted lines: central 
regions are to the right, outskirts are to the left). As long as 
a few per cent only of the gas have been turned into stars, a 
track evolves vertically as the gas surface density does not de- 
crease significantly. Following its vertical leg, however, the track 
bends leftwards, thereby embodying the combined stellar mass 
increase and gas mass decrease. We note that a similar plot is 
provided in fig. 13 of IGutermuth et all (1201 ll) . They assume a 
star formation law where the star formation rate per unit area 
shows a power-law dependence on the gas column density. It 
thus differs from our model which builds on the time-evolution 
of volume densities. 



2.5. Distribution of YSO surface densities 



iBressert et al.l d2010l) recently obtained the distribution of lo- 
cal surface densities of Spitze r-detected YSOs in star-forming 
regions of the Solar Neighbourhood (distance smaller than 
500 pc). Their observed distribution is well-approximated by a 
Gaussian function with a peak at 22 YSOs pc" 2 and a standard 
deviation of 0.85 in log 10 2,YSOs (see their fig. 1). A fully con- 
sistent comparison with the result of IBressert et al.l (1201 Ol) still 
requires our model to be extended to entire molecular clouds 
(i.e. to cumulate many molecular clumps) since their observed 
distribution encompasses a dozen of local star-forming regions. 
It is nevertheless interesting to see what stellar surface density 
distribution our model predicts at its current stage of develop- 
ment. 

Building on the stellar surface density profiles derived in 
Section 12.41 Fig. [4] shows the distribution of the logarithmic 
stellar surface densities, that is, the projected mass of stars en- 
closed within an annulus of radius s, dm p +" 1 {t,s), as a function 
of the annulus logarithmic surface density, log 10 (E*(f, s)). Units 
are M© and M B ■ pc~ 2 . The mass, radius and star formation ef- 
ficiency per free-fall time of the molecular clump are as previ- 
ously. The lines with open symbols are the model predictions for 
the clump density indices po and times t given in the key. In the 
bottom panel, two models are depicted for po = 1.7: t = 1 Myr 
and t = 2 Myr (green curves with open pentagons and circles, 
respectively). The low surface-density regime of each model is 
fitted by a Gaussian (black dotted/dash-dotted lines), with mean 
log 10 £* = 0.55 (f = 1 Myr) and log 10 £* = 0.85 (f = 2 Myr), and a 
standard deviation of 0.71. They illustrate the growth of the stel- 
lar mass and stellar surface density as time goes by. All models 
are bell-shaped, with wider distributions for steeper density pro- 
files, as expected. This suggests that the shape of the observed 
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distribution could be used as a probe into the density profile of 
star cluster parent clumps. 

The Gaussian with which lBressert et all (l2010h describe their 
data is shown as the bottom panel solid (red) line. Assuming 
a mean stellar mass of O.5M per YSO, the surface density at 
their peak O 22 YSOs pc~ 2 ) equates with log 10 Z* = 1.04. Their 
standard deviation is 0.85 and our normalization is arbitrary. At 
t = 2 Myr, our po = 1.7 model is in good agreement with the 
distribution observed in the Solar Neighbourhood. In particu- 
lar, the observed distribution agrees with our model prediction 
better than with the surface density dis tribution of sink parti- 
cles of the hydrodyna mics simulations of Bonnell. Clark & Batd 
(12008b (see fig. 2 in iKruiissen etatl |2012|) . It remains to be 
seen h ow including heavily crowded regions in the lBressert et alJ 
d2010t) 's sample would affect the comparison. Access to high- 
surface density regions (e.g. the core of the Orion Nebula 
Cluster) - which Spitzer fails to resolve and which there- 
fore affects the high-density tail of the observed distribution 
- would actually allow us to test our model more extensively. 
Finally, we emphasize that the observation of a smooth dis- 
tribution in log 10 E* alone does not allow one to conclude 
that star formation is not made of m ultiple discrete modes 
dPfalzner. Kaczmarek & Olczak 1 12012b. For a set of low-mass 
gas-poor clusters, iGieles. Moeckel & Clarke ] J20T2l) propose 
that the surface density at the peak of the observed distribution 
is driven by the degree of early cluster expansion. 

In this section, we have developed a star-forming region 
model building on a constant star formation efficiency per free- 
fall time in a spherical molecular clump with a radial volume 
density profile po oc r~ Po . When po- 1 .7 as observed, the model 
predicts a local star-formation law oc with a - 2. In the 
next section, we combine the model to a surface density thresh- 
old imposed by the stellar background against which the star- 
forming region is seen projected. We demonstrate that this leads 
to the time s equence for emb edded-cluster development recently 
proposed by iPfalznerl (1201 ll) . 



3. The observed growth sequence of embedded 
clusters 

At first glance, the model presented in Se ction [2] departs from 
the scenario proposed by IPfalznerl (201 1 ). In th e scenario she 

' lLada & Lada (120031) . an embedded 



devised based on the data of 
cluster grows outwardly, with outer shells of stars added with 
time. That is, star formation is delayed in the outer regions com- 
pared to the inner ones. In contrast, in the present model, star for- 
mation proceeds all through the molecular clump albeit at a rate 
slower in the outskirts than in the c entre. As quoted in SectionQ"! 
however, the relation inferred by IPfalznerl (1201 ll) between the 
radius and the mean volume density of embedded clusters, i.e. 
Peel K r~ 1 , is reminiscent of one of constant mean surface den- 
sity, i. e. p ec i oc r~~\. This suggests that the data of lLada & Ladal 
d2003l) are surface-density limited. The refore, prior to c ompar- 
ing our results to the growth sequence of lPfalznerl (1201 ll). we in- 
clude in our model this observational bias. Note that I Allen et al.l 
(120071) made a similar point. 

The top and bottom panels of Fig. [5] show the rising with 
time of the volume- and surface-density profiles of the stellar 
component for the model of Section [2] (Mo = 10 4 M o , R = 6pc, 
€ff = 0.1) with a density index po = 1 . 7. The bottom panel als o 
includes the surface density cut-off of ICarpenter et al.l (|2000): 
%bck = 40 M ■ pc~ 2 (horizontal dotted line). The shaded area 
visualizes the surface density regime with < l^ck- We note 



that the embedded clusters of lLada & Ladal (120031) are spread 
over distances from the Sun ranging from =s lOOpc to 2.4 kpc. 
Each cluster has therefore its own limiting background sur- 
face density. Since these are unknown, we resort to the single 
value of ~Lbck quoted above. It should be considered as a mean 
value used for illustrative purposes. The surface density thresh- 
old Hbck = 40 Mr ■ pc~ 2 f ollows closely the embedded-cluster 
track defined by Pfalzner (l2011l her fig. 2). We will come back 
to this point in Fig. [6] 

The surface density limit imposed by the stellar background 
conceals the 'wings' of the stellar component, whose observed 
limiting radius is therefore smaller than the actual one, that 
is, Sbck < R- These observed radii, Sbck> are depicted as (red) 
diamonds in the bottom panel of Fig. [5] They defin e an out- 
ward time-sequence, in agreement with lPfalznerl d201 ll) . In other 
words, while the whole stellar component follows a 'vertical' 
evolution, i.e. the density increases at fixed radius, the fraction of 
it denser than the stellar background mimics a 'horizontal' evo- 
lution, i.e. the observed limiting radius increases with time. In 
what follows, we refer the part of the stellar component 'above' 
the background as the 'embedded-cluster' (non-shaded area in 
bottom panel of Fig. [5] > 'Zbck when s < Sbck)- That is, 
Sbck = fed with r ec i the embedded-cluster radius. We define the 
embedded-cluster mass, m ec i, as the stellar mass enclosed within 
the three-dimensional radius r ec [. 

For each time-step, we obtain the radius and mass of the em- 
bedded cluster and its mean volume density p ec i = 3 • m ec i /(4 • n ■ 
r 3 ; ). Figure [6] provides three different perspectives of the evo- 
lution of the whole stellar component (red open circles) and of 
its associated embedded cluster (blue open circles): from top to 
bottom, volum e density agains t radius (parameter space identi- 
cal to fig. 2 in IPfalznerl 1201 ll) . mass against radius, and mass 
against time since star-formation onset. The size of the open cir- 
cles scales with the logarithm of the stellar mass. 

In the top panel, the vertical (red) sequence at a radius 
R = 6 pc depicts the growth with time of the mean volume den- 
sity of the whole stellar component. Times in Myr and stellar 
masses in M are given to the left and right of the sequence, 
respectively. The decreasing (blue) sequence on the left is the 
predicted embedded-cluster sequence, i.e. p ec i versus r ec i, with 
times since star-formation onset and embedded-cluster masses 
quoted below and above it. This sequence runs parallel to the 
surface density threshold Ebck (solid red line) applied to the star- 
forming region. The difference in volume density between the 
threshold I.j, c k and the predicted sequence (r ec i, p ec i) stems from 
the density profile of the stellar component, that is, the density at 
the limit r ec i = s^ck is lower than the mean volume density inside 
r ec \. The (green) dashe d line depicts th e observed embedded- 
cluster track defined by iPfalznerl (1201 ll) with which our model 
(blue open circles) agrees reasonably well. In a follow-up paper, 
we will map the parameter space (Mo, R, po and e//) to investi- 
gate what patterns emerge in the t~eci~Peci space as a function of 
the input parameters. 

The middle panel of Fig.|6]illustrates the stellar mass growth 
of the whole star-forming region and of the embedded cluster 
(red circles at constant radius to the right and blue circles to 
the left, respectively). Model time-spans t are indicated along 
the sequences. The bottom panel shows the stellar mass in de- 
pendence of the time t since star-formation onset. As the stellar 
mass builds up, an ever greater stellar mass fraction makes it 
above the background limit (see also bottom panel of Fig. [5}. 
Therefore, the shift between the embedded-cluster and stellar- 
component sequences tightens. While the embedded cluster rep- 
resents about one-tenth of the stellar component total mass at a 
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time t — 0.1 Myr, the cluster mass fraction has risen up to about 
one-third by t = 2.5 Myr. We therefore conclude that, not only 
does the stellar mass fraction in clusters depend on the surface 
density threshold adopted to define clusters (as opposed to their 
surrounding 'haloes' of stars), it also depends on when the star- 
formation process started. 

Here, a remark about the embedded-cluster mass function 
is worth being made. The mass function of star clusters is a 
powerful diagnostic tool of their evolution. For instance, if the 
mass function slope remains unchanged between the onset and 
the end of violent relaxation 0, this implies that cluster infant 
weight-loss through vio lent relaxation is mass-ind ependent (e.g. 
iParmentier et al.L [20081). The embedded clusters of lLada & Ladal 



(120031) have a power-law mass function of slope =a -2 (their 
fig. 2), similar t o the mass function of young gas-free clusters 
(IChandar et al.L 12011 . We now see that this aspect could be 
co nsidered under a new viewing angle. If the clusters compiled 
by lLada & Lada (120031) mark different stages of their forma- 
tion process, as iPfalznerl j20Tlb suggests, then the slope of the 
embedded-cluster mass function at gas expulsion (i.e. when the 
build-up of the stellar content terminates) may differ from -2, 
since the mass of the clusters at the start of the growth sequence 
should be corrected for the still missing stellar mass. This would 
have consequences as to whether cluster infant weight-loss is 
mass-independent or not. 

For an assumed volume density profile of the stellar com- 
ponent, one can extrapolate the mass of the embedded cluster to 
that of the whole stellar component by representing the star mass 
distribution in the r - m(r) space. Figure [7] illustrates the rela- 
tions between the distance r from the star-forming region centre 
and the enclosed mass m(r), for the gas initially and the stellar 
component at t = 2 Myr (symbol/colour-coding identical to pre- 
viously). For a pure power-law density profile of index po, radii 
and masses follow (taking the case of the initial gas density dis- 
tribution, Eq.Q}: 

mp(r) _ I r \ 3 "Po 
Mo \RI 

Note that a steeper volume density profile shows up as a shal- 
lower r-m{r) relation (e.g. in Fig. [7] the gas tracks are steeper 
than their star analogs). The dash-dotted (black) straightline has 
a slope of 0.5, as expected for the stellar component in the low- 
density regime when po = 1.7, i.e. 3-q^3-3po/2- 0.5. While 
it indeed fits well the model with po — 1 .7 in the clump outer re- 
gions (open circles at r > 1 pc), the deviation increases towards 
the clump centre as the stellar density profile gets shallower than 
q = 3po/2 (see middle panel of Fig. [TJ- The lines with open 
symbols provide a direct mapping of how much stellar mass 
is enclosed within a given radius. For instance, when po = 1.7 
and t = 2 Myr, the embedded cluster has a radius r ec i ^ 1 pc and 
a mass m ec i =s 400 M , while the stellar mass enclosed within 
fl^6pc is > 10 3 M o . 

Finally, we note that the time t of our model corresponds to 
the age-spread of the model stellar component. It differs from 
the mean age of the stars which would be inferred by an ob- 
server. Figure [8] illustrates the mass-weighted mean age of the 



5 Violent relaxation is the dynamical response of the stellar com- 
ponent to the expulsion of the unprocessed star-forming gas. Its dura- 
tion - from =! 1 Myr to several tens of Myr - depends on the crossing- 
time of the cluster gaseous precursor and on the observed cluster re- 
gion. The denser the gas clump (i.e. the shorter its mean crossing- 
time), the smaller the aperture with which the cluster is observed (i.e. 
the higher the observed cluster de nsity), the faster the cluster evolution 
dParmentier & BaumgardtTl2012h 



stars against the time-span t of the simulations. In case of a con- 
stant star formation rate, the mean age is simply 0.5f. However, 
in our model, star formation slows down with time and the first 
stars to form (i.e. the oldest ones) tend to dominate the mass 
budget. The mean age is therefore higher than 0.5f, although the 
deviation is moderate, at most 10% of half the simulation time- 
span. The effect is stronger for po =1.9 than for po =1.5 since 
a steeper gas density gradient slows down star formation more 
than a shallow one (see Section B.l.ll and Fig.|9]below). 

The top x-axis of Fig. [8] shows the time t since star forma- 
tion onset in units of the initial local free-fall time at the clump 
half-mass radius, r/, m , when po = 1.7. Using Eq. [2T| we find 
r/ lm =i Q.6R for po = 1.7, which leads to an initial local density 
Po(fhm) - 12M Q ■ pc~ 3 (Eq.[TJ and an initial local free-fall time 
Tf/(t = 0, rhm) — 2.4 Myr. A scaling in units of the initial free-fall 
time enables us to apply our model to molecular clumps with dif- 
ferent mean densities and, therefore, different rates of star forma- 
tion (since the mass of newly formed stars depends on the ratio 
t/Tff{t = 0,r); see e.g. Eq. [5J. For instance, a molecular clump 
that is 100 times denser initially forms stars at a rate that is 10- 
times faster (i.e. the bottom x-axis of Fig. [8] shrinks by a factor 
of 10). This density-dependent star formation rate is akin to the 
density-dependent rate of dynamical e volution of star clusters af- 
ter gas expulsion studied in depth by IParmentier & Baumgardfl 
(120121) . 

4. Model Consequences 

4.1. Star Formation Efficiencies: Global and Local 

4.1.1. Global Star Formation Efficiency 

Figure [9] shows the evolution with time of the global star forma- 
tion efficiency, SFE, namely, the stellar mass fraction averaged 
over the whole molecular clump: S FE - M+/Mo. The density 
indices of the precursor clump are identical to those in Fig. [1] 
i.e. po - 1.5,1.7 and 1.9 (see key). The symbol-free lines de- 
pict the tangents to the models at the onset of star formation. 
A steeper density profile concentrates a higher gas mass frac- 
tion into the clump inner regions, which have higher densities 
and shorter free-fall times. This accelerates star formation ini- 
tially (i.e. the tangent for po - 1.9 is steeper than for po = 1.5) 
and leads to higher global star fo rmation efficiencies (see also 
iTan. Krumholz & Matznerl 120061) . With the ongoing gas deple- 
tion, star formation slows down as the comparison between the 
symbol-free lines and the numerical models illustrates. We note 
that, by a time t = 2.5 Myr, the global SFE remains low, of or- 
der 0.1-0.15. We will discuss this result in the framework of star 
cluster survival in Section l4~2l We also remind that the predicted 
SFEs are parameter-dependent, here obtained with Mo = 1O 4 M , 
R = 6pc and e// = 0.1. Higher clump masses or smaller radii 
would lead to higher star formation efficiencies at a given time t. 

4.1.2. Radially-Varying Local Star Formation Efficiencies: 
Two- and Three-Dimensional 

To define the three-dimensional local star formation efficiency 
at the distance r from the clump centre, we build on the local 
volume densities of stars and gas: 



_ Mr) 
P*(r)+Pg(r) po(r) 



(22) 



We remind the reader that the subscript '0' refers to the gas 
clump prior to star formation. Equation [22] is illustrated in the 
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top panel of Fig.[10]for the same model as previously at a time 
t — 2 Myr (symbol/colour-coding identical to Fig. [9]). As already 
introduced through the top panel of Fig. Q] the local star for- 
mation efficiency is higher in the clump centre than in its out- 
skirts. It is worth keeping in mind that the derived efficiencies de- 
pend on the adopted star formation efficiency per free-fall time, 
on the assumed duration of the star-formation process, and on 
the clump mean volume density. That is, a higher e/f, longer 
time t and/or a larger clump mean volume density 2>MqI(AtiR?) 
would all increase 63o(r). Also note that steeper density profiles 
achieve higher values of e^o since the clump central density is 
then higher. 

f3o(r) is not the star formation efficiency inferred by ob- 
servers, however, since observers work with surface densi- 
ties rather than volume densities. We therefore define a two- 
dimensional star formation efficiency based on the observed sur- 
face densities: 



£2d(s) = 



(23) 



Equationl23lis shown in the bottom panel of Fig.fTUl 

At a given radius s — r, eio(s) is lower than 63£>(r). This is 
so because the star surface densities in the clump centre vicinity 
intercept outskirt material where the initial gas volume densities 
and achieved star formation efficiencies are lower. This is further 
illustrated in Fig. Q~T] which depicts eio in dependence of 630. 
The dotted (black) line corresponds to eio = Qd, while the solid 
one obeys eio - -0.20. The solid line thus shows that, for the 
case of relevance here, the difference between the - measured - 
two-dimensional and - actual - three-dimensional star formation 
efficiencies can be as high as 20 per cent. 

Both 63 £> and e 2 D rise from a few per cent at the clump edge 
up to higher than 50 per cent at the clump centre. In particu- 
lar, star formation depletes almost the entirety of the initial gas 
content in the clump centre and reaches values as high as 60- 
90 % in the inner r < 0. 1 pc. This contrasts with the low global 
star formation efficiency, SFE, which we saw in Fig. [9] This 
has important consequences for the survivability of (part of) the 
stellar component as a bound cluster after residual gas expulsion, 
as we shall discuss in Section [4721 The high central star forma- 
tion efficiencies also allow us to propose an alternative explana- 
tion as to the low gas content in the central r egions of the most 
ev olved molecular clum ps of IHiguchi et al.1 (120091) . According 
to IHiguchi et al. l d2009T) . this absence or scarcity of gas results 
from its ongoing dispersal. In contrast, in our picture, the ab- 
sence of gas in the central regions of deve loped embedded clus- 
ters (e.g. the Orion Nebula Cluster) or of IHiguchi et al.l d2009l) 
clumps stems from most of it having been fed to star formation, 
even before gas dispersal starts. 

The top panel of Fig. [12] shows how the three-dimensional 
star formation efficiency increases from a few per cent to almost 
unity as the initial gas volume density rises from a few M Q ■ pc -3 
(near the clump edge) to higher than 1O 4 M • pc~ 3 (toward the 
clump centre). We stress again that other model parameters will 
be conducive to different efficiencies. With po — 7M - pc -3 = 
«H 2 100 cm -3 andpo >7000 M -pc~ 3 Efi H , > 10 5 crrT 3 , where 
«h 2 is the molecular hydrogen number density, the range of vol- 
ume densities characterizing the molecular clump under scrutiny 
here extends from the diffuse mole cular gas probed in C 12 Q up 
to the dense gas traced i n CS (e.g. Shirley et all 120031) or dust 
continuum emission (e.g. lSchuller et aUl2009l) . 

The bottom panel of Fig. Q~2] illustrates the 2-dimensional 
star formation efficiency in dependence of the residual gas 



surface density, that is, e 2 D versus E g . This is the relation 
which would be inferred from observing the star-forming 
region. T he two filled squares highlight the measurements 
made by iGutermuth et al.l d201ll) : e 1D = 2.3 % and 26 % at 
Eg = 20 and 300 M e ■ pc~ 2 , respectively. They are practically 
consistent with our model despite it not being based on detailed 
star and gas density profiles. It is interesting to note that 
star formation efficiencies of order 1% happen at gas surface 
densities only slightly higher than the limit beyond which 
the gas becomes entirely molecular in spiral galaxies, i.e . 
2 = 9M Q • pc~ 2 (see top left panel in fig. 8 of iBigiel et all 
120081) . Therefore, the lowest star formation efficiencies would 
correspond to the lowest possi ble surface densities where the ga s 
is entirely molecular. [But see lLada. Lombardi & Alves 1 d2010b l 

Finally, Fig. [T3]is the analog of Fig. |9]for the local star for- 
mation efficiency, of which it shows the evolution with time for 
different radii r and density indices po. While the star formation 
rate is nearly constant in the clump outskirts (i.e. the local star 
formation efficiency is linearly increasing), the severe gas deple- 
tion in the clump inner regions slows down star formation there. 

4.2. Radially-Varying Local Star Formation Efficiency and 
Cluster Survival 

Figure [9] shows that at the end of our simulations, the star for- 
mation efficiency averaged over the whole star-forming region 
remains low, SFE =* 10-15%. At first glance, this is signifi- 
cantly smaller than the threshold required for the stellar com- 
ponent to retain a bound star cluster after residual gas expul- 
sion. The relation between SFE and the fraction of stars remain- 
in g bound to the stellar compon ent, Fbound, is shown in fig. 1 
of iParmentier & Gilmord (120071) based on the sim ulations of 
iGever & Burkertl d200ll) and iBaumgardt & Kroupal d2007l) . For 
the sake of clarity, this Fb oun( ]-vs-SFE relation is reproduced as 
the solid (black) line in our Fig. Q~5] It assumes instantaneous 
gas expulsion and a weak external tidal field. The threshold to 
retain a bound group of stars (Fbound > 0) is SFE t i, 0.33. We 
stress that this threshold necessarily depends on the hypothe- 
ses of the models used to derive it. The vast majority of cluster- 
gas-expulsion models developed so far assume that the local star 
formation efficiency is uniform all through the star-forming re- 
gion, that is, 63d is independent of r and constant (as a result 
SFE=63d). This is in stark contrast with our result of Section l4~Tl 
that the local star formation efficiency spans near to two orders 
of magnitude, from almost unity down to a few per cent (top 
panel of Fig. ITOb . This property of the star-forming region has 
important implications for the ability of the stellar component 
to retain a bound star cluster after gas expulsion. Owing to its 
small residual gas fraction, the clump central region has a high 
resilience to gas expulsion. It will therefore necessarily produce 
a bound cluster. Most of the gas-expulsion-driven disruption of 
the stellar component takes place in the outskirts where the local 
star formation efficiency is low. 

The top panel of Fig.[T4lshows the local star formation effi- 
ciency in dependence of the enclosed stellar mass fraction. The 
bottom panel depicts e^o against the normalized radius r/R (i.e. 
another representation of the top panel of Fig. \W\ with the x- 
axis scale now linear). Let us consider the model with a den- 
sity index po = 1.9 (blue line with open triangles). About 15% 
of the stellar mass arises from gas having experienced a local 
star formation efficiency 630 > 0.60. Under the assumption of a 
uniform star formation efficiency in that limited central region 
(0.60 < 63D < 0.95 when r/R < 0.04 or r < 0.25 pc; see bottom 
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panel of Fig. IT4t. we can use the solid black line in Fig.[]3]to in- 
fer the bound fraction of stars after violent relaxation. More than 
85% of the stellar mass remains bound in the central region. We 
can now derive a lower limit on the mass of a surviving bound 
cluster. With a global SFE of 13% at a time t = 2 Myr (Fig. 0, 
the whole stellar component has a mass M* ^ 1 300 M , of which 
15% are located in the central area practically unaffected by gas 
expulsion. Therefore, a cluster of at least 200 M e will remain 
after violent relaxation in spite of the low global SFE. In other 
words, over the global scale of the whole star-forming region, 
the bound fraction F bound must be at least 15%. 

lAdamsl d2000i) developed a pioneering model accounting for 
a radially-dependent local star formation efficiency, which quan- 
tifies the effect introduced above. That is, because of the high 
local star formation efficiency in the central regions (i.e. high 
€}d at small r), a cluster survives gas expulsion (Fbound > 0) 
despite a low global star formation efficiency (low SFE). The 
relations he derived between SFE and Fbound are shown in 
Fig. Q3] as the dash-dotted (red) and dashed (blue) lines. They 
correspond to isotropic and anisotropic stellar velocity distri- 
butions, respectively. For SFE =a 0.1-0.15, uniform-e3£> mod- 
els predict the stellar component full disruption (solid black 
line: F b ound - 0). In contrast, the radially-varying models of 
lAdamsl d2000t) give bound fractions as high as 20-50%. With 
Mo = 10 4 M o , SFE= 0.13 and Fboumi - 0.3, our star-forming re- 
gion would thus give rise to a bound star cluster with a mass: 
Ma - Fb oun d -SFE ■ Mo - 400 M . This is in fair agreement with 
the lower limit of 200 M we derived above. 

5. Discussion and Future Work 

In this section, we firstly discuss how our results and those 
from hydrodynamics simulations of star formation relate to each 
other. We then touch on two aspects of our work which we will 
refine in the near future, namely, the value of the star formation 
efficiency per free-fall time (e//), and the impact of the density 
profile of molecular clumps. 

5.1. Clump density profile and embedded-cluster 
morphology 

iMtiller et alj J2002b detected in their sample of star-forming 
molecular clumps a correlation between the clump radial 
density profile and mass, with higher-mass clumps being on 
the average steeper (top panels of their fig. 18). The observed 
p-index range (0.7 < p < 2.5) can now be related to the 
comprehensiv e set of adaptat i ve-me sh refinement simulations 
performed by iGirichidis et alj d201 lb . These simulations show 
that the initial gas density gradient of spherical molecular 
clumps is of paramount importance to the morphology of the 
star clusters they form. In clumps initially uniform in density, 
supersonic turbulence has the time to compress the gas in locally 
disconnected areas before global collapse sets in. This results in 
locally disconnected filaments and spatially distinct subclusters 
of sink particles. In contrast, centrally concentrated clumps (e.g. 
a Bonnor-Ebert or power-law po = 1.5 density profile) form 
filaments more strongly connected and centrally concentrated 
because of the shorter time-scale for global collapse in the 
clump central region. This eventually f avours the form a tion o f 
one main central cluster (see fig. 4 in IGirichidis et al.L 1201 ll) . 
The morphology of these 'single-block' clusters - substructured 
or smooth power-law de nsity profile - d e pends o n the initial 
turbulent velocity field (IGirichidis et al.L 1201 2bl) . Solenoidal 
turbulence in Bonner-Ebert or po = 1.5 density profiles can 



produce embedded clusters with a smooth radial density profile. 
For instance, the embedded-cl uster Q-paramet e r H ca n be 
of order 1.2 (figs 11 and 12 in IGirichidis et all 1201 2bl) . We 
note that 0.9 < Q < 1.5 corresponds to a power-law density 
index 2 < q < 2.9 for the emb edded cluster (see Table 1 in 
ICartwright & WhitworthL |2004|) . a g-range which agrees well 
with our model predictions (see our Eq. [8] q < 3/?o/2 = 2.25 
if po = 1.5). We caution, however, that our model does not 
account for the clump global collapse and are therefore not fully 
comparable to those of Girichidis et alJ d201 ll) . Finally, we note 
that these adaptative-mesh refinement simulat ions predict that 
the firs t stars form in the clump central region (IGirichidis et all 
l2012al) . Once secondary stars form around the central objects, 
they accrete the infalling gas preventing it from reaching the 
clump central region ('fragmentation-induced starvation'). This 
forces new protostars to form at increasingly large dista nce from 
the clu mp centre, an aspect which agrees well with iPfalznerl 
(2011) 's scenario. 

If the correlation be t ween the clump density index and mass 
found by iMuller et alj d2002l) is genuine (a bias arising from 
more massive sources being located at larger distances where 
the resolution is poorer cannot be fully excluded yet), one would 
expect shallow low-mass clumps to produce ensembles of sub- 
clusters, while steeper high-mass clumps would give rise to 
'single-block' clusters. Which mode of star cluster formation - 
subclustered or 'single-block' - dominates the process of star 
formation then depends on the slope of the clump mass func- 
tion. O bserved mass functio n slopes are shallower than -2 (e.g. 
- — 1 .7, iRramer et al. I fl998h . High-mass clumps therefore domi- 
nate the star-forming gas mass and most newly formed stars may 
thus originate from 'single-block' clusters. Low-mass clumps 
and subcluster ensembles would dominate the cluster population 
only in terms of number. 

Sim ilarly to the uniform density profile of IGirichidis et alj 
(1201 lh. the end product of t he hydrodynamics simulations of 
i Bonnell. Clark & Bate! d2008l) is an ensemble of subclusters (see 
llVIaschberger etalj, |2010L for their merging history). This is not 
surprising since the gas cylinder they scrutinize is character- 
ized by a weak density gradient along its main axis. Despite 
a low global SFE (SFE 0.15), their subclusters are gas-poor 
on a length scale of 0.1 -0.2 pc, an effect due to the accretion 
of gas ont o sink particles and th e accretion-induced subcluster 
shrinkage dKruiissen et al. L I20L21) . This is reminiscent of the lo- 
cally high SFE achieved by the central region of our molecular 
clump despite a globally low one (Figsl9land[T0ll. Note, however, 
that the respecti ve embedded-cluster morpholo gies are different 
(subclustered in IBonnell. Clark & Bate] d2008l) vs. a smooth ra- 
dial density profile in our model). Given their high local SFE 
and small number of stars, the relaxation time of subclusters is 
short enough for their early dynamical evolution to be dominated 
by collis ional stellar dynamics rather than violent relaxation. As 
a result, iMoeckel et alj d2012l) find that subclusters expand due 
to the scattering of t heir stars, eventual ly erasing the subclus- 
ter structure (see also lSmith et al.Ll2008l) . We remind that these 
results hold for the low-mass end of the cluster mass spectrum 



6 The g-parameter ( ICartwright & Wmtworthl 120041) quantifies the 
degree of substructures in star clusters. It is defined as the ratio be- 
tween the normalised mean separation of the stars and the normalised 
mean length of the edges of the minimal spanning tree associated to 
the (proto)cluster. A cluster with Q < 0.8 is substructured/fractal, while 
Q > 0.8 indicates a smooth radial density profile 



10 



G. Parmentier and S. Pfalzner: Local-Density Driven Clustered Star Formation 



(i.e. cluster mass of the order of 100M Q ), typical of the Solar 
Neighbourhood. 

Finally, we emphasize that the outcomes of smooth-particle 
hydrodynamics and adaptative mesh refinement simulations of 
star-forming regions necessarily depend on their input physics. 
It remains ill-known how the inclusion of magnetic fields or ra- 
diative feedback (hence gas heating), how treating the proto- 
stars as extended gas spheres rather than sink particles, would 
affect the emerging properties of model star-forming regions. 
Additionally, these simulations do not currently access the 
regime of ab initio formation of high-mass clusters. This is up to 
at least two orders of magnitude higher than the gas and stellar 
masses covered so far. How relevant to the formation of high- 
mass clusters the physical processes at work inside low-mass 
clouds are remains an open question. 

5.2. The star formation efficiency per free-fall time 

To match the observ ed local star formation law of 
iGutermuth et all (1201 lh with our model (M - 10 4 M o , 
R =s 6pc) at t — 2Myr requires a star formation efficiency per 
free-fall time of e/y = 0.1. We insist that the so-derived value of 
€ff is parameter-dependent. A longer time t or a higher clump 
density (i.e. shorter free-fall time) would both lead to higher 
YSO surface and volume densities. That would require a smaller 
e/f to maintain the good match between the model and the 
scaling law averaging the observations shown in the top panel 
of Fig. ft 

At this stage, it is therefore premature to draw any conclu- 
sion about our e/f e stimate being an order of magnitude higher 
than th at derived by Ikrumholz & Tanl d2007l) . IKrumholz & Tan! 
(120071) consider the fraction of the Milky Way star formation 
rate taking place in different classes of objects, from Giant 
Molecular Clouds (number densities «h, - 100cm -3 ) to HCN- 
traced molecular clumps («h 2 - 6 ■ 10 cm" 3 ). Building on the to- 
tal mass of these objects in the Galaxy and on their mean density, 
they derive a star formation efficiency per free-fall time of about 
one per cent, independent of volume density (SFR/y 0.01 in 
their notation; see their eq. 1 and fig. 5) . On the oth e r hand , in 
a study of individual molecular clouds, lEvans etal I J2009h in- 
fer 0.03 < SFR/y < 0.0 6 (their section 4.3). This is intermediate 
between the estimate o f IKrumholz & Tanl (120071) and ours. 

Ideally, to derive e/f in the framework of our model, we 
need the total masses of gas and stars enclosed within a given 
radius (not simply orders of magnitude as used in this intro- 
ductory paper), as well as the time elapsed since the onset of 
star formation. In a forthcoming paper, we will apply our model 
to smaller-scale molecular clumps with known size, gas mass 
and star mass. This will allow us to put on a firm footing 
the co mpari son between our e tf esti mate, those of lEvans et all 
(120091) and IKrumholz & Tanl (12007b . The comparison will be 
especially interesting as these three works define a sequence 
of distinct spatial scales, from Galaxy-integrated star formation 
(IKrumholz & Tanl 12007 ). to star formation in individual molec- 
ular clouds ( lEvans et al 



20091) and star formation in individual 



molecular clumps (our work). 



5.3. The time-evolving gas density profile of molecular 
clumps 

Summing up the observed surface density profiles of the stars 
and (unprocesssed) gas in star-forming regions will provide a 
direct estimate of po, the initial density index. In that respect, 



our assumption of a single pp-v alue valid over 6 pc in radius 
is probably an oversimplification. iPirogovl (120091) finds that the 
density index of massive-star-forming clumps is p — 1.6 ±0.3 
within 0.8 pc from th eir centre. Beyond tha t distance, the density 
profile drops steeper. iBeuther et all (|2002|) reach a similar result 
for the massive star-forming regions they study. Two power-laws 
are needed to describe their clump density profiles: p — 1 .6 + 0.5 
within 32"from the clump centre, and a steeper density profile 
beyond (their fig. 2). The break at 32"corresponds to half of their 
finable range and equates with 1 pc at the typical 6 kpc distance 
of their star-forming regions. Therefore, it will be interesting in 
the future to convert our present one-zone model into a two-zone 
one and to quantify how this affects the estimate of e/f. 

In ad dition, it is worth noting that the densit y in dices re- 
ported bv lBeuther et al.l(l2002l) . lMuller et al.l(l2002l) and lPirogovl 
(120091) refer to star-forming molecular clumps. Yet, the bottom 
panel of Fig.Q~]shows that, as the stellar content builds up, the gas 
density profile gets shallower, again reflecting the faster gas-to- 
star conversion in the clump inner regions. Over the radial range 
0.2-1.0pc and over a time-span of 2.5 Myr, this gas density pro- 
file evolves from po = 1.7 to p ^ 1.3. It is therefore likely that 
observed star-forming molecular clumps had at star formation 
onset steeper profiles than is observed now. Models starting with 
po = 1.9 may thus be more appropriate than those with po = 1.5 
or po - 1.7. One way of testing this observationally would be to 
compare the density pro files of starless molecular clumps (e.g. 
iTackenberg et al.Ll2012l) and star-forming ones. 



6. Summary and Conclusions 

We have presented a model quantifying the stellar content of 
molecular clumps as a function of time and initial gas volume 
density. The model key-ingredient is the star formation effi- 
ciency per free-fall time, e/f, namely, the mass fraction of gas 
turned into stars per free-fall time, r// (Eq. |2}. The model orig- 
inality resides in it building on a local free-fall time defined in 
relation to the local volume density of the star-forming gas. That 
is, the radial volume density gradient of molecular clumps leads 
to a radially-dependent free-fall time, shorter in the clump cen- 
tral regions and longer in the clump outskirts. In other words, 
star formation proceeds more quickly at shorter distance from 
the clump centre. As a result of this differential rate of star for- 
mation, the radial density profile of the stellar population built 
by the clump (which we refer to as the 'stellar component') is 
steeper than that of the clump initially (i.e. q < 3po/2, where q 
and po are, respectively, the density indices of the stellar com- 
ponent and of the clump initially; see Eq.[8]l. 

We combine the volume density profile of spherical molec- 
ular clumps (Eq. [TJ with their radially-dependent free-fall time 
(Eq. |2]i to model the density profile of the stellar component as 
a function of time and gas initial density. This is done both nu- 
merically (Section 1272b and analytically (Secfion l2~3l >. 

We then project the volume density profiles for the unpro- 
cessed gas and stars and obtain the corresponding surface den- 
sity profiles (Section 12.4b . Under the assumptions that YSOs 
do not migrate away from their birth sites significantly, and 
that molecular clumps do not experience significant gas mo- 
tions and/or inflows and/or outflows, predicted surface density 
profiles are compared to the observations. Spe cifically, we con- 
sider t he local star formation law inferred by IGutermuth et all 
(1201 lh . It relates the local surface density of molecular gas, "L g , 
and the local surface density of YSOs, Z*: X* =a 10~ 3 Sg, with 
a =* 2 and the surface densities in units of M • pc' 1 . We find 
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that a stellar density profile steeper than a gas density profile, 
as predicted by our model (Eq. [7] and top panel of Fig. [0, nat - 
urally leads to a > 1, as observed by iGutermuth et all (1201 ll) . 
In particular, molecular clu mp density indices pq in t he range 
1 .5 < po < 1.9 dBeuther et al.Ll2002tlMuller et alll2002l) are con- 
ducive to a =s 2, in excellent agreement with what is observed 
(top panel of Fig. [3]). 

To compare the normalizations of the observed and predicted 
star formation laws, we consider the specific case of a molecu- 
lar clump with an initial gas mass Mo 10 4 M o and a radius 
R =* 6 pc. These mass and radius are chosen so as to emulate the 
most prominent concentration of YSOs in the MonR2 molecular 
cloud (Section l2T2b . We find that the predicted star formation law 
agrees with its observed counterpart about 2 Myr after the onset 
of star formation (i.e. t = 2 Myr) for a star formation efficiency 
per free-fall time e/f =0.1. For that same model, we derive the 
YSO surface-density distribution (Section [2.51 Fig. |4j) and find 
a good agreement with th e observed dis t ributio n of the Solar 
Neighbourhood inferred bv lBressert et alJ (l2010h . We stress that 
this comparison is only preliminary given that our present model 
covers one single m olecular clump, while the distribution of 
iBressert et alJ d2010h encompasses several star-forming regions. 

In the second part of the paper, we combine our model of a 
time-evolving star-forming region with a surface density thresh- 
old (Section [3]). Prior to the one-to-one identification of YSOs 
based on their infrared excess with the Spitzer telescope, the 
subtraction of a given surface density accounting for the con- 
taminating foreground and background stars was often relied on 
to define the limiting radiu s of a cluster and its mass (see sec- 
tion 2 of lAllen et al. 1 120071 for a summary of the different tech- 
niques for cluster member identification). The combination of 
our model predictions with a surface density threshold can thus 
be compared with pie-Spitz er cluster data sets, s uch as the star- 
clu ster catalog com piled bv lLada & Ladal (120031). 

iPfalznerl d201ll) finds that the clusters of lLada & Ladal (12003b 
define a sequence of decreasing volume density along with in- 
creasing radius. She suggests this sequence to be a time sequence 
corresponding to the growth of the cluster stellar content, an hy- 
pothesis we confirm. As time goes by, the total stellar mass in- 
crease raises the fraction of the stellar component seen above 
the surface density threshold (we refer to this 'emerging' part 
as the 'embedded cluster'). Therefore, the observed embedded- 
cluster radius gets larger too (red diamonds in the bottom panel 
of Fig.0. The mean surface density of the embedded cluster re- 
mains approximately constant since the cluster definition rests 
on a surface density limit. A constant surface density equating 
with p ec i oc f~\, the observed growth of the embedded-cluster 
radius thus leads to a decrease of the observed mean volume 
density. In contrast, the outer radius of the stellar component as 
a whole stays constant and its mean volume density increases 
with time (see the sequences of red and blue circles in the top 
panel of Fig. [6j. 

In Section |U we have discussed the impact of a radially- 
varying local star formation efficiency on the post-gas-expulsion 
evolution of clusters. That the density profile of the stellar com- 
ponent built by the clump is steeper than that of the gas initially 
implies a local star formation efficiency increasing towards the 
clump centre. As a result, the clump central region is more re- 
silient to gas expulsion than the clump outskirts. Despite a low 
global star formation efficiency (SFE 10-15% at f = 2 Myr for 
Mq =2 10 4 M o , R^6pc and e/y =0.1, see Fig. [9]), the stellar com- 
ponent can still leave a bound cluster after violent relaxation due 
to the high local star f ormation efficie ncy fen) in the clump in- 
nerst regions (Fig.fTOt. Adams! d2000h computed the bound frac- 



tion of stars at the end of violent relaxation for a centrally -peaked 
star formation efficiency profile. He found that such a config- 
uration increases the cluster survivability compared to models 
where the star formation efficiency is uniform all through the 
molecular clump (Fig.[T"5l). 

Although our model addresses the formation of leaky (i.e. 
low-density) clusters, it already provides a hint as to why 
some starburst (i.e. high-density) clusters may be in virial 
equilibrium despite their very yo ung age (e.g. Westerlund- 1 , 
iMengel & Tacconi-GarmanL 120071) . The top panel of Fig. [T2] 
shows that volume densities characteristic of starburst clusters 
(i.e. several 10 3 M o -pc" 3 at least) lead to three-dimensional 
star formation efficiencies in excess of 50% (for the star for- 
mation efficiency per free-fall time adopted here). A high star 
formation efficiency limits the impact of gas expulsion on the 
cluster dynamics and the departure from virial equilibrium af- 
ter gas expulsion. [Besides, a high density also implies a short 
crossing -time which hastens the return of the cluster to equi- 
librium dParmentier & BaumgardtLl2012l) l. That does not imply 
that gas expulsion is not a major driver of young cluster dynam- 
ics in general. Only high-density regions - e.g. dense molecular 
clumps close to the Galactic centre and central regions of mas- 
sive molecular clumps in the Galactic disc - may be relatively 
immune to gas expulsion. 

Finally, we note that a sound quantification of the stellar 
component also requests an accurate estimate of the star for- 
mation efficiency per free-fall time (e//) since this parameter 
drives the relation between the local star formation efficiency 
63D and the local volume density po(r) at a given time (top panel 
of Fig. [121). Should e// be halved, the predicted global and local 
star formation efficiencies will be reduced by a factor of two (at 
most). Here arises a direct connection between the rate of star 
formation in molecular clumps and the stellar dynamics of post- 
gas-expulsion clusters. That is, the rate of dissolution of young 
clusters and the star formation efficiency per free-fall time are 
tightly related. 
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dius R = 6pc, density index po = 1.7, a star formation efficiency 
per free-fall time e/f = 0.1 and the times t quoted in the key. 
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Fig. 14. Top panel: Relations between the local three- 
dimensional star formation efficiency , ejoif), and the enclosed 
stellar mass fraction, m*(< r)/M+. Bottom panel: £3o(r) against 
normalized radius, r/R, where R is the outer radius of the star- 
forming region. Symbol/colour-coding identical to Figs.l9lfT3l 
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